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Abstract 

We investigate the quantum origin of the primordial cosmological gravitational waves 
in the Brans-Dicke theory in the two conformally related frames, the Jordan-Brans-Dicke 
frame and the Einstein frame. We calculated the theoretical observable in both frames 
and we compared both with General Relativity. We compute the number of gravitons Nk 
produced during inflation and the observables: power spectrum Pj, spectral index ut and 
energy density Q-k. The comparison shows that for the case of the particles number Nk 
the results are the same in both frames and in General Relativity when the Brans-Dicke 
parameter is much bigger than unity. For the spectral index nj we show that it is possible 
to get a scale invariant perturbation in the Jordan-Brans-Dicke frame when CO oo and in 
the Einstein frame when ft) ±0°. In both frames, the results found for the power spectrum 
Pj and the energy density Q. show that the preferred values of CO are differente from that 
that are found in the local tests. 

KEYWORDS: scalar-tensor gravity, gravitational waves, conformal transformations 

PACS numbers: 04.30.-w, 98.80.-K 

1 Introduction 

The gravitational waves are tensorial fluctuations in the metric of spacetime. This particular 
perturbation is not explicitly coupled with the energy density and pressure of the matter of the 
Universe and does not contribute to the gravitational instability that generates the cosmological 
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structures we see today. On the other hand this study is of great interest because it supplies the 
specific signature of the metric theory of gravity. There is a special interest in the primordial 
gravitational waves of quantum origin because they may have left special signatures in the 
polarization of the Cosmic Microwave Background (CMB) anisotropics [lj-;3J. Moreover, since 
the gravitational waves were generated in the early Universe, before the time of last scattering 
of CMB photons, they can be a window to the primordial phase of evolution of the Universe and 
can help testing the initial conditions of the scalar field. These waves are predicted by Einstein's 
theory of General Relativity (GR), but they still have to be directly detected. Great efforts have 
been done in this sense and there is hope that a new generation of experiments in space may 
allow this detection ffl. The gravitational waves spectrum frequency extends over a wide range 
of interest, from 10^ //z to lO^Hz, depending of the sources that generate those waves. Many 
works have been done in order to identify specific signatures of cosmological models in the 
spectra of gravitational waves, for instance, in the case of quintessence model [5] and string 
cosmology [^-^8J. 

In 1999, the SN la observations [j9| showed that the Universe is currently undergoing accel- 
erated expansion. A possible theoretical explanation for this acceleration is the vacuum energy 
with negative pressure, called dark energy (that violates the strong energy condition). There are 
many dark energy models, for example quintessence [ |10J , k-essence [11 1, phantom [ |12| and 
Chaplygin gas [J3J. The simplest candidate for the dark energy is the cosmological constant A 
whose equation of state is p — —p. But there are some problems with this alternative as, e.g., 
why the cosmological constant is so small, nonzero, and comparable to the critical density at 
the present? So, it is interesting to study dark energy models as a possible theory of k-essence 
coupled to gravity since it involves the simplest form of nonlinear kinetic term for the scalar 
field (j) [fT4). 

The existence of a classical scalar field in the nature has been considered in many alternative 
theories of gravitation. Scalar-tensor gravity theories are a particular case of the huge class 
called Extended Theories of Gravity which are revealing a useful paradigm to deal with several 
problems in cosmology and related areas including fundamental physics [ |15J . The prototype 
of this kind of theories is the Brans-Dicke theory (BDT) ||T6- 19 1 and its revival leads us to ask 
if the generation of gravity waves can be modified through the introduction of this scalar field. 
Here we propose to study the evolution of tensorial fluctuation in the traditional BDT and its 
quantum origin and calculate the observables that can be compared with the future observational 
data. The main interest in the presence of the scalar field is the possibility that it can change 
the evolution of gravity waves in comparison with the results obtained with GR [20). But the 
situation is not so simple: there are other problems related with the BDT as the reference frame 
problem. 

As others scalar-tensor theories, the BDT can be formulated in the Einstein frame (EE) or 
in the Jordan-Brans-Dicke frame (JBDE), which are conformally related. We have two con- 
formally related frames, ds^ and ds^, if we multiply the line element ds-^ by a nonvanishing 
function of spacetime coordinates Q?-{x^) and the resulting line element have the same light 
cone that the first one. Thus, the purpose of the conformal transformation is to give another rep- 
resentation of the theory that is equivalent to it. This can be made by a simple transformation of 
units [pTJ. In the case of the BDT we are looking for a system of units where the gravitational 
coupling is constant. 
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In the specific case of the BDT which conformal frame is the physical one is a very con- 
tentious issue p2| . The most common argument against the JBDF found in the literature is that 
the scalar field energy density can assume negative values. This happens because of the non- 
mininal coupling between the scalar field and geometry. The terms with the second covariant 
derivatives of the scalar field in the Brans-Dicke field equations contain the connection and 
therefore part of the dynamical description of gravity. But the energy density can be made non- 
negative with the help of a new connection f JJ^, [j23| which is also associated with the physical 
metric g^v 

By a conformal rescaling of the metric g^v = g^vigjiv^^) and a nonlinear field redefinition 
^ = ^(g^v, ^) we obtain the EF of the BDT, that has a positive definite energy. In this frame, 
however, there is a violation of the equivalence principle, due to an anomalous coupling of the 
scalar field ^ to ordinary matter. This violation is small and compatible with the available tests 
of the equivalence principle. It is, indeed, a low- energy manifestation of general compactified 
theories where, e.g., the dilaton in the n-dimensional Kaluza-Klein theory plays the role of the 
Brans-Dicke scalar field ^. 

Classically both the theories are equivalent and describe the same physical system with 
different names. The equivalence appears as follows: the space of solutions in the JBDF is 
isomorphic to that the EF and the isomorphism is given by the relations g^y gfiv and ^ — )• ^. 
But from the physical point of view the situation is very different mainly due the in definition 
concerning the energy and related to the coupling between the scalar field and the gravitation. 
Moreover, another factor that distinguishes the two frames from the classical view is that in 
the case of EF gravitational radiation is generated by the quadrupole term while in the case of 
JBDF we also find the terms of monopole and dipole radiation [24|. Nevertheless it has been 
a tradition to consider the Jordan metric as physical in the BDT, since the initial proposal of 
Brans-Dicke. On the other hand, a fundamental assumption in the theory of gravitation is that 
the gravitational interaction is generated by the massless spin-two graviton. Once one accepts 
this hypothesis, one must accept the EF as the physical spacetime metric in the BDT. Therefore, 
so as we can see, this is an unsolved problem related to the BDT. To shed light on this problem, 
far from theoretical issues, we calculate the physical observables related to the BDT through the 
quantization of gravitational waves and we hope that the observations decide what is the best 
answer to this important question. 

This Letter is organized as follows. In section 2 we study the background solutions, the 
gravitational waves equation, the quantum gravitational waves equation and the physical ob- 
servables in the JBDF. In section 3 we make the same analysis but now in the EF. Finally, in 
section 4 we discuss the results and show our conclusions. 

We use in this Letter standard tensor notation: the greek indices run from zero to three; 
the latin indices run from one to three; the spacetime metric is taken to have signature +2, or 
(— , -f-, -|-, -|-); covariant derivatives are indicated by semicolons; partial derivatives are indicated 
by commas and the scalar field ^ is a time function. Moreover, we work using units in which 
h — c — kB — I. 
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2 Jordan-Brans-Dicke frame 



2.1 Background solutions 

In this scenario, the action of the BDT is given by p5| 



5^= d^x 



(1) 



We assume the background of the Universe as being spatially flat (k = 0), homogeneous and 
isotropic, i.e., it is described by the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

ds^ = -dt^ + a^{t){dx^ + dy^ + dz^) , (2) 

where a{t) is the scale factor of the Universe. This choice is because the present observations 
indicate that the density parameter Q. is close to unity p6| . With the metric in this form we are 
considering a state on a certain three-dimensional space-like surface in spacetime, which is due 
to the topological structure of the FLRW spacetime - the usual canonical ADM decomposition 
of the metric. 

Considering the matter content of the Universe as a barotropic fluid, p = ap, we obtain the 
equations of motion 

~'-a = JP[ 3^ j+^^n ■ 

The above equations must be supplemented by a conservation equation of the energy- 
momentum tensor 

p+3%{l + a)=0 . (5) 

Background solutions can be obtained assuming that the scale factor a{t) and the scalar field 
(t) have a power-law form 

a(f)oc?'- , (l){t)<^f . (6) 

In what follows we shall consider the inflation case (a = — 1). In this scenario we have a 
drastic expansion of the Universe during the early period of the Big Bang. The background 
solutions are 

s = 2 , r=a + - . (7) 



2.2 The gravitational waves equation 

Cosmological gravitational waves are obtained by means of a small correction hij in equation 
(|2]), which represents the metric. Hence, the general expression of the metric ([2]), related to the 
unperturbed metric, is replaced by 

d^2 = _cl?2 + [a^ (t) 5ij + hij] dx'dxj . (8) 



4 



We use the above metric in the field equation to obtain the resulting gravitational waves equation 



2 2— 



h,j = 



(9) 



where k is the wavenumber, the dots indicate cosmic time derivatives and we have written 
hijit.x) = h(t)Qij, where Qij are the eigenmodes of the Laplacian operator, such that Qu = 
Qki,k = 0. 

Performing the transformation from the cosmic time to the conformal time, a{t) dr\ = dt, 
and representing the derivatives with respect to r\ by primes, equation (|9]) assumes the form 



r- 2— + 2^-2^ 



= 



a 

To resolve this differential equation we need to perform the Fourier transformation 



(10) 



(11) 



where the polarization tensor e^-^ (^) can be decomposed as e^-^- (^)e'-'(^')(p) = b'^^ik — 
p) to the two polarization states e ©. Then, we have 



a 2 (j) 4\^y 



M(A)(^7)=0 



(12) 



In this way, we rewrote the equation ( 10) in terms of a harmonic oscillator 



We can express the equation (|l2fin terms of the new parameter 4>(t]) = a{ri) ^y (j) [r]) 

,2 ^"(ri) 



Hence, the equation ( 12 1 becomes 



4.(7]) J 



M(A)(^7) =0 



(13) 



With the background solution of the scale factor and the scalar field and considering a = — 1 
we get 

4>"(t]) 2(3 + 2(0) (1+2(0) 

-^:r^ = ^. — ^ \o I 19 . with -oo<T]<-T]i . (14 
4>(t]) (1 -2(0)2 |t]|2 ( - (1 

The solution in the inflationary era is 



li{kr]) = y^{Aij4'\k\7l\)+A2j4^\k\7^\) 



(15) 



where Ai and A2 are integration constants, Jif^^^ and J^^^^ are the Hankel functions of first and 
second kind, respectively, and V = 2{i^2(o) order of the Hankel functions. 
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2.3 The quantum gravitational waves equation 



Since the primordial gravitational waves are of quantum origin we need to quantize the per- 
turbation hij. In order to do that we need to put the equations for the perturbation into the 
Hamiltonian form. Dirac p7| showed that it is possible to make the quantization without the 
problems resulting from the diffeomorphism invariance. In order to proceed with the quantiza- 
tion we start with the Lagrangian density of the gravitational waves in the BDT, given by 



167t 



(16) 



that leads to two corresponding Hamiltonians 

1 



-2 / 



4)'' 



11^ + d'jiijiij 



(17) 
(18) 



where n and ft are the canonical momenta. The above Hamiltonians differ each other by addic- 
tion or subtraction of a total conformal time derivative. 



The expression of the quantum Hamiltonian obtained by the equation (18) for the gravita- 
tional waves in the BDT is constructed by using the canonical commutation relations 



[fl{x,ri),7t{y,ri)] = id'^^\x-y) 



(19) 



So, by transforming the normal modes of oscillation of the gravitational waves in field 
and momentum operators in the Heisenberg description of quantum mechanics we obtain the 
Hamiltonian 

4)" 



(2) 



d^x 



(20) 



where the field operator is constructed from its classical field as follows 

hij -^hij{x,ri] 



1 f£lf\k)d^k. 



(2;r)3/2 



with the operator /i^ given by 

and the function satisfying the following equation 



(21) 



fk+ k 



4) 



fk = 



(22) 



(23) 



which is the same differential equation that describes the behavior of the classical function 
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Now, to obtain the coherent states of the relic gravitons in the BDT we consider the creation 
and annihilation operators 



(24) 



whose temporal evolution is given by the Heisenberg equation and by the Hamiltonian (17), 
with the general solutions given by 



a%{ri) = Uk{ri)a-^{riQ)+Vk{ri)aK{rio 



where rjo is a fixed initial time and the functions Uk{r}) and Vyt(T/) behave as 



(25) 



-— = ikuk + —Vi, 
dT] " ^ k 



-— = ikvk + —ui. 



(26) 



The equations (25) are the Bogoliubov transformations [28 1 of the gravitational waves in 
the BDT and Uk{ri) and Vk{r]) are the Bogoliubov coefficients that satisfy the relation 



\Vk\ 



1 , 



(27) 



result that guarantees the unity of the temporal evolution of these operators. 

The number operator is defined as = v^v^. So, in our model the number of gravitons is 
given by 



Nk = \vk\ 



\fk 



2k i^^nl^ 



+ 



1 

2k 



dfk 



dr] 



l^/^,dJ7 dA\ 1 
2k^y'dr]'^'^'dnJ 2 



(28) 



where the relation between and is 

* i dfk , fk 



When kr]» I the number operator simplifies to 



k — i 
2k \ ^ 



1 /4> 



2k^ V ^ 



1 

2 



(29) 



(30) 



If A;T] << 1 we have 



1 



(31) 



In Figures [T]j2] we can see the behavior of the number of particles in the inflation era in the 
BDT in the case of the JBDF to different values of ft) and small and big values of the frequency 
and we compare with the EF and with the GR. 
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2.4 Physical Observables 



In this section we apply the formulation of quantum gravitational waves obtained in BDT in 
the JBDF in order to determine the most common and interesting observables related to the 
gravitational waves. The observational data in the present do not allow us to include or exclude 
the presence of a primordial spectrum of relic gravitons compatible with any model and any 
frame. The idea is that these theoretical results calculated here might in the future be compared 
with the observations obtained by the projects related to the search for gravitational waves of 
cosmological origin. 



2.4.1 Power Spectrum 

The power spectrum is related with the two-point correlation function of the tensor modes in 
the following way p9| 



{0\hUx,rj)h'Jiy,r^)\0)= / 77)^^d(lnfc) , 



kr 



(32) 



where Prik.r}) is the power spectrum and the state |0) is annihilated by a? 



This operator is given by the equation pT] ) where the vacuum state |0), in which there are 

(33) 



no particles, is annihilated by the annihilation operator 

a_i\0)=0 . 

The calculations of the expectation value of the field operators yield 

d'k \fkm\-ik.r.- 



{0\hl{x,r])h'j{y,r])\0)^64K 



(2;r)3 



y) 



(34) 



Knowing that gravitons can be produced in isotropic models and using spherical coordi- 
nates, we have 



{0\hl{x,r])h'jiy,r^)\0)=647r 



k" |Ap2sin(fcr) 
[InY kr 



dk . 



(35) 



Comparing the equation ( [32] ) with the previous equation we have the power spectrum of the 
tensor modes, within the JBDF, parametrized as 



PT{k,r\) =?>2k 



l_\fil_ 
71^ a^^ 



(36) 



In the inflation phase the function f]^{k\r]\) is defined as 



fkim) 



r]\^i^\k\7)\) 



and by using the small argument limit of the Hankel functions |30|, we have 

8;r22vr2(v) k^ 



Prik^t 



{m\) 



l-2v 



(37) 



(38) 
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By definition the Hubble parameter is given by H 
this model can be written as 

'10} + l\ 1 



So, the conformal time in 



(39) 



2(0-1) aH 

Thus, the power spectrum in the BDT during inflation for the modes out of the horizon becomes 



PT{k,r]) 



8;r22^r2(v) f2co + l 



2(0-1 



1-2V ^2 



\aH 



3-2v 



(40) 



Li Figures [3j|7] we show the behavior of the power spectrum Pj {k, 7] ) as a function of the 
conformal time |t] |, the frequency /(Hz) and of the parameter (O. 



2.4.2 Spectral Index 

The spectral index is obtained by the expression 

dln/'r(fc,7?: 



(41) 



where a = aH = k is the moment when the tensorial modes cross the horizon again. For 

this calculation we introduce the parameter e, that measures the rate of decrease of the Hubble 
parameter during the inflation phase and it is given by 



1 



1 dH 



In the JBDF we obtain that 



Thus, we have 



riT = 



dt \ H J aH^dri ' 
2 



2(0 + 1 



dink 



In 



= 2 



k dH d\n(j)\ 
'H~dk~ dink J 



(42) 



(43) 



(44) 



The first term on the right side is calculated as follows 



^kdH 
^H~dk 



k dH d|T]| 



Hd\ri\ dk 



k dH fd\ri 



and, by using the equation (39), we have 



d|77| 



dk 



d_ 
die 



2(0 +l\ 1 



2(0-1 J aH 



H dr] \ dT] 



-i- 

dk\k 



-1 



d\ri\ 



dk 



(45) 



1 

1? 



(46) 
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In this way, with the help of the equations ( [43| ), ( [45] ) and ( [461 ), '^he equation ( [44] ) can be 
written as 

'2co+l\ dln0 



-2e 



2(0-1 



din A: 



(47) 



When (0 ^ oo and considering ^ as a constant we obtain the same result that in GR, /. e. , 



tiT = -2e [31 1 



Likewise, we can calculate the term dependent on the scalar field (j) in equation ( [47] ) 

1 



dink 



kd(j) 




k d^ dr] 




k d0 


fd|77l 


^ dk 


a 


^ dri dk 


a 


" ^dri 


[dii 



d\ri\ 



dk 



(48) 



that with the help of equation (39) and with the background solutions of the scalar field we 
obtain the following result for the spectral index nr 



riT = -2e 



2(0 + 1 
2(0-1 



2(0-1 



(49) 



Again, the above result gives the same value of the spectral index of the tensorial modes of GR 
at the time of inflation, when (O — t- oo. 

We can write the expression of nj in terms of the parameter (O. In order to do that we replace 



the relation ( [43] ) in the result found above, such that 

-8 



2(0-1 



By using (42) we obtain 



2(0 = - - 1 , 
£ 

which can be replaced in the final result for the spectral index such that 

4e . _i 



£-1 



-4£(l-e)" 



(50) 



(51) 



(52) 



For small values of e we can expand the term ( 1 — e) ^ obtaining 

nr = -4e(l+e + ^(£2)) , (53) 
and considering only zero-order terms we have 

tiT ^ -4e . (54) 

2.4.3 Energy Density 

From the Lagrangian density ( [16] ) we obtain the energy-momentum tensor of the gravitational 
waves 



(55) 
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Considering 



we find that 



where 



A=igi,© 



(56) 



(57) 



(58) 



The energy density is the time component of the energy-momentum tensor of the gravita- 
tional waves 

p = TQ''^(^(^+g'^dihdjh] . (59) 



To calculate the quantum energy density we first consider the following redefinition 

h-^h= — j= , 
ay/0 

such that the expectation value of the energy density is given by 

(p) = (0| P |0) = ^ ((0| h'h'* |0) + (0| • V^* |0)) . 



(60) 



(61) 



With help of the expansions 



(2;r)3/2 



(2;r)3/2 



^{i^)e-''-y+7q{i^)e''-y 



the solutions 



A^ = fl^(T?o)A(T?)+a^(T]o)A*(n) , 



(62) 



"k 

and the commutation relations 



%'4 



(63) 



(64) 



we obtain the expectation value of the energy density of the gravitational waves in the BDT that 
can be write as 



(P) 



1 r k^dk 



a^^J (2;r)3 



2 

lAI --^{gkfk+gkfk) 



2sin/rr 
kr 



(65) 
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where we have considered the presence of an isotropic background of relic graviton and r 
\x — y\. In the limit ^ — )• y we have 



(P) 



a 



1 fk^dk 



1%^ 



(66) 



The spectral energy density per logarithmic interval of the wavenumber is defined as 

1 d(p) 



pc dlnfc 



(67) 



I. e.. 



Sn^-^ k^ 
3Hl iTt^a"^ 



Ifkl^-^igkfk+fkgl) 



(68) 



where pc is the critical density of the Universe. 

For all the oscillation modes that are inside the Hubble radius (z. e. kr] » 1), we have 
^ ^ikfi^. This implies that 



1 /4)'V 



and, rewriting it in terms of the power spectrum we have 



(69) 



^l{k,r^ 



k^ ^ (, 1 

;Pr 1 + 



2k^ \ 4) 



For modes outside the Hubble radius (z. e. fcr/ << 1) we have gk 
us to find that 



Snk^ 



■]fk\ 



which can also be rewritten as a function of the power spectrum Pj 

k^ 



(70) 

^fk- This result allow 
(71) 



(72) 



We stress that when the modes re-enter the horizon the energy density is proportional to the 
power spectrum. The behavior of the spectral energy density in the JBDF is shown in Figure [8] 



3 Einstein frame 

In this section we calculate de general behavior of the BDT but now in the EF. As in all scalar- 
tensor theories, the action of the BDT is preserved under a group of field redefinitions that 
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contains two functional degrees of freedom. Specifically, if one defines a new metric g^y and a 
new scalar field ^ by 



(73) 



where = e^^ is a function of the coordinate and C7 is a new scalar field. We also consider 
the following definitions: 

2 



K 



2(0 + 3 



(74) 



3.1 Background solutions 



The conformal transformation (73 1 change the theory's frame from the JBDF to the EF and the 
action can be written as 



d\v^ 



R 1 



(75) 



where this Lagrangian was obtained with the help of an integration that is necessary to eliminate 
a boundary contribution [32|. To obtain the fields equations we have to variate the action (175]) 
with respect to the metric tensor and the scalar field o, and we get 



1 



1 



□ C7 



(76) 
(77) 



where fuv is the stress-energy tensor in the EF. It is related with the stress-energy tensor in the 
JBDF T^v by the expression 7)iv = ^^^'^^T^v- 
The stress-energy tensor conservation gives us 



(78) 



Applying the fiat Friedmann-Robertson- Walker metric in the field equations (76) and in the 
conservation relation ([78]) and considering a perfect fluid with equation of state p = YP getting 



y 
1 



b 

C7 + 3-<7 = -Ae 
b 2 

p+3^(i + r)p 



'^"(37-l)P 

3 . 
: -A(7(l + 7)p 



(79) 
(80) 
(81) 



where b{t) is the scale factor in this scenario. The equation (81 1 can be easily integrated to give 
us 



pb 



3(i+r)g-i(i+r)(^ 



(82) 
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where C2 is an integration constant. 

In order to solve the equations ( [79] ), ( [80] ) and ( [8T] ) we will assume that the field a and the 
scale factor b have the following power law solutions: 



For the case of inflation, when 7 = —1, the solutions are 

1 



1 , P± 



1±A 1 + 



16/3^ 



P=C2 



(83) 



(84) 



In the following we consider the p+ solution because it reproduces the inflationary model in 
the limit —°°<T}<—rii. Using the conformal time transformation dt = b{ri)dr} , where t and 
T] are the cosmic and the conformal times, respectively, we have 



b{ri)=bo[bo{p+-\){-ri)y-p+ , -oo<T]<-T]i , 
where — r/i is the time when the radiation era begins. 



(85) 



3.2 The quantum gravitational waves equation 

We can perturb the metric tensor such that = V^v + h^v, where f]^v is the background 
metric tensor and h^v is the metric perturbation. Using the perturberd metric tensor in the field 
equation (|76|) we obtain the gravitational waves equation in the EF: 



b- /V b , - 

h,,--h,,+ (--^-2-]h,, = , (86) 



where V = ri'^didj. 

In order to solve the above equation it is convenient to make the transformation 

where k is the wave number and the polarization tensor e^j^ (k) obeys the relation ej-j^ (k)e'-'^^"> ik') = 
25 (k — k')d^^''''^ . The wave equation is then 



where 11' = dji/drj and 



with 



b" F(j82 



M=0 , (88) 



b rj 



2 



(89) 
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To solve the wave equation in this frame we make the definitions /i = y^g{z) where z = kri 
to obtain 



and the solution is 



^{kr]) = ^{Aij4'\k\ri\)+A2J^^r'{k\ri\)) , 



(2)/ 



(91) 



(92) 



for V = v/^(/32) - 1/4. 

Following the quantization method described in [3J we transform the field jl into operator 
of the Heisenberg description of Quantum Mechanics: 



The coefficients fi^ can be written as 



(2;r)3/2 



(93) 



(94) 



where T]o is a fixed comoving time and the functions fkif]) and the operators d-^{Tjo) obeys the 
relations 



fkfk ~ fkf'k — ' ' 



(95) 

(96) 
(97) 



We have already found the solution for equations like equation (95). To have a well behaved 



function in the limit r\ — )■ —oo and with help of the normalization relation (96) we find that for 
inflation 



fk 



(98) 



where 0^ is a phase factor. Still following the prescription of [[3| the number of gravitons is 
given by 



Nk 



Ifkl 
2k 



k' + 



f\ 2^ 



+ 



1 

2k 



dfk 



dri 



2k b y'dr] ^^'drj 



1 

2 ' 



(99) 



In Figures [T]j2] we show the behavior of the number of particles in the inflation era in the 
BDT in the case of the EF to different values of O) and small and big values of the frequency 
and we compare with the JBDF and with the GR. 
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3.3 Physical Observables 
3.3.1 Power Spectrum 

The quantization process gives us 

the same result obtained in GR. So, the gravitational waves power spectrum is given by 

|2 

The behavior of the power spectrum is shown in Figures |3|jT 



4^2:^^ . (101) 



3.3.2 Spectral Index 

In order to find the spectral index of the gravitational waves during inflation in the EF we have 
first to calculate the parameter e of inflation, given by 

d^ 1 d^ 1 

where H = b/b is the Hubble parameter in the EF During this calculation process, with help of 
the relation bH = ^ ^ it is possible to find 



that is necessary to calculate the spectral index. 
The spectral index is given by 



dlni'r(^,T?) 

riT = 



dink 



-2e-^ = — , (104) 

1-P+ 1-P+ 



where a = bH = (^p^j k and the power spectrum Pt is evaluated for krj << 1. We found 
that when (o ±oo the spectral index goes to 2e differently from the JBDF where, in the same 
limits, ht — 2e. 

3.3.3 Energy Density 

The energy-momentum tensor of the gravitational waves is 
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Using the quantization ( |100[ ) we find 

1 rk^dk' 



(P) = (0|p|0) 



Ifkl^-jiSkfk+glfk) 



(106) 



where gj^ = and |0) is the vacuum state. 

The spectral energy density per logarithmic interval of the wavenumber is defined as 



Q{k,7i] 



1 d(p) 
pc dink 



I. e.. 



k^ 



\gkV+\k' + 



b' 



\fk\ 



b' 



(gkfk+fkgl) 



where pc is the critical density of the Universe. 

The density energy behavior is shown in the EF in Figure |9j 

4 Conclusions 



(107) 



(108) 
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Figure 1 : The number of gravitons produced by inflation in the BDT in the different frames and 
in the GR as a function of the frequency /(Hz) with fix conformal time |t/ | = 10^^ for different 
values of CO, including negative values in the EF. 
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f = 10'^ Hz 




2.x 10"^ 4.x 10"' 6.x 10"^ 8.x 10"' 1.x 10"'' 



1^1 

Figure 2: The number of gravitons produced by inflation in the BDT in the different frames and 
in the GR as a function of the conformal time |t] | with fix frequency / = 10^^ Hz for different 
values of CO, including negative values in the EE 

In this work we considered the quantum gravitational waves production during inflation in 
the BDT in the conformally related JBDE and EE We saw that, in contrast with the GR, not 
only the scale factor but also the Brans-Dicke scalar field models the time evolution of the 
cosmological gravitational waves in both frames. In the JBDE the scalar field appears in the 
gravitational waves equation and in the EE the scalar field changes the time evolution of the 
background. 

We computed the quantum particles number Nk created during inflation and we saw that 
the number predicted by this model in the JBDE approximates that predicted by GR when the 
Brans-Dicke parameter co is big enough (depending on the frequency, as we can see in Eigures[T]- 
[2]). The approximation between the EE and GR is achieved when the parameter co is negative 
and big. 

Afterwards we calculate three observable parameters: the quantum power spectrum Pj^, the 
quantum spectral index nj and the quantum energy density T]) of the gravitational waves 
in both frames. The result obtained for the power spectrum i\ shows that, different from the 
local observations and from the previous result acquired for the graviton number A^^, this model 
in the JBDE approximates the GR result one if CO assumes values not so big, like ft) = 10, at 
least for high frequency. Eor low frequency higher values of the Brans-Dicke are preffered, like 
CO = 10^, as can be seen in Eigure[3j Eor co > 10^ the curve moves away from the GR one for all 
frequencies. In the EE we can see (Eigure|4]) that this model approximates GR for small values 
of CO, specially for high frequency. 

Eor the spectral index we obtain a different result. We can compare the outcome of the GR 
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Figure 3: The power spectrum of the gravitational waves in the BDT and GR, in logarithmic 
scale, as a function of the frequency /(Hz) with fix conformal time \ti\ = 10~^ and different 
values of ft) in both frames. 
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Figure 4: The power spectrum of the gravitational waves in the BDT and GR as a function of 
the frequency /(Hz) with fix conformal time |t] | = 10~^ and small values of ft) in both frames. 
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Figure 5: The power spectrum of the gravitational waves in the BDT and GR as a function of 
the conformal time \r] \ with fix frequency / = 10^^ Hz and large values of CO in both frames. 
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Figure 6: The power spectrum of the gravitational waves in the BDT and GR as a function of 
the conformal time |r] | with fix frequency / = 10^^ Hz and some values of (O in both frames. 
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f=10^^ Hz 




Figure 7: The power spectrum of the gravitational waves in the BDT and GR as a function of 
the conformal time \r}\ with fix frequency / = 10^^ Hz and some values of (O in both frames, 
including negative values in EE 
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Eigure 8: The graviton energy density in the BDT in the JBDE and in GR as a function of the 
frequency /(Hz), for small values of (D and fixed = 10~^. 
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Figure 9: The graviton energy density in the BDT in the EF and in GR as a function of the 
frequency /(Hz), for negative values of CO and fixed jr/ 1 = 10^^. 



for the spectral index nr with that of the BDT in the JBDF when CO >> I, result that agrees 
with the one obtained for the quantum particle number A^^, and when (O — )• oo it yields a scale 
invariant perturbation. In the EF we can get scale invariant perturbation too, but for both limits 

CO — > ±°o. 

In the case of the quantum energy density big values of CO make the energy density to 
explode and the GR prediction is approximated when co is of the order of unity in the JBDF and 
in the EF ft) < but not CO < —20 is preffered, depending on the values of the frequency and the 
conformal time. 

The results obtained here suggest that the coupling constant CO of the BDT can vary with 



the cosmic scale [33,34| and it depends substantially on the chosen frame. When gravitational 
waves are detected we can put some constraints on CO for large scales. 

Until there we can compute the scalar perturbations during inflation using the BDT, compare 
the results with these from the GR and, with the available observational data, calculate de 
ratio rsT between the tensor and scalar power spectra. The ratio rsr is related with the power 
spectrum of the polarization modes of the CMB and, combined with the spectral index of the 
scalar perturbation, it allows us to rule in or out our inflationary model [35| 36 1. We can also 
perform these same calculations but using from the beginning a background metric that is not 
flat but with a curvature, i.e., k^O. Other possibility is to do the same for both scalar and tensor 
perturbation but in a model where the Brans-Dicke parameter co varies with the scalar field 0, /. 
e., CO = C0{^). These calculations are in current investigations and will be presented elsewhere. 
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